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Abstract 

We present a study for a design of a long baseline appearance ex- 

periment to probe the high sin^ 1Q and low Am^ region relevant to explain 
the atmospheric neutrino anomaly. The experiment relies on a good iden- 
tification of quasi-elastic interactions, which is a clean topology that has 
an important contribution in the lowest Am^ part of the region probed. 
The detector we studied is a fine grained liquid scintillator detector of 
15 kilotons, optimized to detect electrons from r ez/eZ^r decays, while 
rejecting backgrounds from 7r° in interactions and electrons from the 
z/e beam contamination. As a reference, the proposed v^^ neutrino beam 
from CERN to Gran Sasso was used. 



1 Introduction 



The question of whether or not neutrinos have a mass has been in the mind of physicists for 
some time. If the neutrinos have non-degenerate masses and are mixed, neutrino oscillations 
could in principle be observed. The study of this phenomenon, by which transitions between 
neutrino flavours are possible, constitutes a very precise way to determine the mass states of 
the neutrinos. 

Although there have been some indications from solar neutrino experiments (Homestake 
GALLEX [0, SAGE [||, Kamiokande and Super-Kamiokande P) that the neutrinos indeed 
could oscillate, the recent results from the Super-Kamiokande atmospheric neutrino exper- 
iment gave the strongest indications, by showing evidence for i/^ disappearance. Other atmo- 
spheric neutrino experiments have obtained results compatible with this hypothesis (Kamiokande 0, 
MACRO § and S0UDAN2 0). Furthermore, the LSND 0] short-baseline appearance ex- 
periment found indications for z/^ oscillating to u^. 

Now that the first indications for neutrino oscillations have been found, the confirmation of 
the oscillation signals and the determination of the oscillation parameters represent important 



objectives. This year already, the K2K |TT| long baseline disappearance experiment will shed 
some light on the situation. At the end of the data taking, it will be able to probe the region 
above Am^ ~ 2 x 10~^ eV^. The MINOS |]12| experiment at Fermilab, which will start taking 



data early next century, also aims to probe the atmospheric region. 



Given the CHOOZ |jT3| result, which excludes the interpretation of the atmospheric neutrino 
data as oscillations above Am^ ~ 10~^ eV^, a plausible interpretation of these data 

is Un Ur oscillations. Proposals for several experiments to detect Ur appearance using the 
proposed NGS fl^ neutrino beam from CERN to Gran Sasso are being prepared. 



Prompted by the recent efforts, both at Fermilab and at CERN, to develop long baseline 
projects with the NuMi |jl2| and NGS beams, the purpose of the present paper is to study the 



possibilities for a simple u^- appearance experiment to probe the Am^ region of the atmospheric 
neutrino experiments. We argue that a totally active scintillator detector in a long baseline 
beam would provide an efficient way to verify the oscillation claim from the Super-Kamiokande 
data. 



The present study is based on the parameters of the proposed NGS beam. This beam has 
a mean energy of about 18 GeV and a baseline of about 735 km, corresponding to the distance 
between CERN and the Gran Sasso laboratory. In order to understand the role of the baseline, 
we also study the performances of the same neutrino detector when positioned 3000 km away 
from the same neutrino source. 

In these conditions, low Am^ values (below 5 x 10~^ eV^) are such that the oscillated Ur 
energy spectrum peaks at low energy (below 10 GeV). At these energies, the quasi-elastic Ur 
interactions have a significant contribution to the total event rate. An experiment aiming to 
reach very low values of Am^ could thus be optimized for the detection of quasi-elastic processes, 
since they represent a very clean signal compared to the deep-inelastic processes. An efficient 
way of studying those interactions is to focus on the events where the produced r decays to 
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> rp, with r — > ez/eZ/,- ). The signature of these interactions is similar to 
u^n ep ) events that 1216 aims to study to search for z/g 



an electron [] ( z/^-n 
the z/e quasi-elastic 

appearance. Figure |I] shows the two topologies involved for z/^ (a) and z/^ (b) 

appearance experiments. The 1216 experiment would be looking for appearance in a z^^ 
beam by comparing the ratio of Vy, and quasi-elastic interactions in detector modules at two 
locations (130 m and 885 m from the target). The current design of the modules consists in 
fine grained fully active scintillator calorimeters amounting to an overall mass of 500 tons. This 
experiment would use the CERN-PS to get a neutrino beam of energy of about 1.5 GeV. In 
the present study, we use one of the detector designs of 1216 (to be described later) and apply 
it to the detection of interactions. 




Vg quasi-elastic 




Vx quasi-elastic 



Figure 1: Schematic comparison of z/gU — > ep and z/^n — > rp with r — ez/gZ^T quasi-elastic 
processes. 



Although addressing the question of neutrino oscillations in the atmospheric region is of 
great importance, the verification of the oscillation claim of LSND is no less important. The 
KARMEN ||T^ experiment is trying to clarify the issue, and other projects have been proposed 
(MiniBoone JTSf and 1216 itself). According to this study, we conclude that the use of a single 
detector technology could allow to search for neutrino oscillations both in the LSND region and 
the atmospheric region by doing a short baseline z/^ — z/g experiment as a first step to a long 
baseline appearance experiment. 

This paper is organized in the following way. A description of the detector is given in 
section |^. The event selection and the backgrounds are discussed in section which is followed 
by section ^ on the oscillation sensitivity. 



2 Description of the detector 



In the studies related to the 1216 Letter of Intent |[I6|, it was realized that a fully active liq- 
uid scintillator detector with a granularity of the order of a few centimetres fulfils the general 
requirements for electron identification and topological reconstruction necessary for a appear- 
ance search. We want to apply the same technique to the detection of z/^n rp interactions 
with subsequent decay r ez/gZ/^, in a detector about thirty times larger (15 ktons). It should 
be noted that liquid scintillator detectors have already been successfully used in a wide range 
of experiments (for example: CHOOZ [|r^ and MACRO 0), and foreseen for future large scale 



applications (for example: KamLAND |[T^ and Borexino [pO|] ) 



^The branching ratio of the r decay to electron is about 18% |T5| 
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2.1 Detector design 



Because of the large scale, a simple and modular structure should be adopted (see figure H). 
The active target consists of an oil based scintillating mixture contained in a large vessel. The 
granularity is provided by optical separators immersed in the liquid. These separators are 
structured in planes of parallel strips. The planes are perpendicular to the beam direction. 
The strips of a given plane are orthogonal to those of the adjacent planes in order to be able to 
do bi-dimensional reconstruction. The light collection is accomplished by wavelength-shifting 
(WLS) fibres placed inside each strip. The light from the fibres is readout by multi-pixel 
photon detectors. The general detector characteristics are reported in table |l] and discussed in 
the following sub-sections. 



Table 1: General features of the detector 



Mass 


15 ktons 


Modules 


4 


Module width 


10 m 


Module length 


42 m 


Sampling 


4 cm 


Granularity 


4 cm 


Planes per module 


1042 


Strips per module 


260417 


Photodetectors per module 


2035 



HPD / PM 



y 




Figure 2: Schematic view of the detector planes 
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2.1.1 The active medium 



Mineral oil scintillator has been studied in great detail, and excellent performances in stability, 
light yield and attenuation lengths have already been reached in the past. Contrary to surface 
readout detectors (for instance KamLAND and Borexino), where the attenuation length of the 
liquid plays an important role because the light must travel a large distance, the light emitted 
by the scintillator in the case of a WLS fibre readout is captured locally inside the fibre. The 
attenuation length of the liquid can then be shorter, which gives more freedom in the choice of 
the mixture. 

A very recent development of oil based liquid scintillator consists in using new solvents 
like PXE (phenyl-o-xylylethane) and LAB (linear alkylbenzenes) Mixtures using 

these solvents can have very good scintillation properties (for example, 87% of anthracene light 
output can be reached with the Bicron scintillator BC599-16 |^). Their physical properties 
are well suited for a very large scintillating WLS fibre detector. They are non-toxic, non- 
flammable (flash point above 145°) and good insulators. Both solvents mentioned here are 
produced in large quantities for industrial applications. These physical properties are crucial 
if the photodetectors are to be placed inside the scintillating mixture without any specific 
precaution concerning high voltages and power dissipation. 



2.1.2 Granularity 

In this study, we fix the longitudinal sampling frequency to 4 cm (corresponding to 1/13 radi- 
ation length) with a transverse granularity of 4 cm, to allow for topological reconstruction of 
the events. With such a granularity, the detector could be placed above ground. According 
to simulation and laboratory tests the proper refiectivity of the walls can be obtained 
in various ways and with different materials (for example: aluminium painted with Ti02, or 
extruded polypropylene). 



2.1.3 Wavelength shifting fibres 



The WLS fibre readout technology allows to collect the scintillating light from a large volume 
(the volume of a scintillator strip) with good efficiency onto a small surface (the fibre diameter 
at one end). This is crucial in our large scale application, because a significant fraction of the 
overall cost is given by the photodetectors. It is then important to maximize the amount of 
active mass which is read by a given photocathode surface. For this reason, we want to have 
strips as long as possible while having enough photoelectrons to measure the ionization loss. 
We assume a strip length of 10 meters and a single-sided readout with a mirror on the opposite 
side of the fibre. From studies done for MINOS |T2[, we expect that a measured signal of an 
average of 10 photoelectrons f\ for a minimum ionizing particle passing through the far end of 
a strip could be achieved with current technology. 



^We assume a typical quantum efRciency for the photodetector (see next paragraph). 
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2.1.4 Photodetectors 



Multi-anode phototubes and hybrid-photodiodes (HPD) are commercially available. They allow 
to read a large number of channels, up to a few hundreds, on a single device. In the case of 
HPD's, the cost is mainly proportional to the surface of the photocathode rather than to the 
number of channels. Typical quantum efficiencies at a wavelength of 520 nm (the peak value 
of the emission spectrum of the WLS fibre) is in the region of 13% to 20%. We have checked 
that these commercial photodetectors can be immersed in the mentioned liquid mixtures. This 
aspect is crucial to minimize the length of the WLS fibres outside the strip and the consequent 
light losses. 

2.1.5 Infrastructure, civil engineering and cost 

The feasibility of a 15 kiloton detector also depends on civil engineering issues like safety, 
assembling, environmental impact and design costs. Industrial oil containers have capacities 
largely exceeding our requests and would be suited for our application. These could then be 
recycled for other purposes after the experiment is finished. The overall detector can be split 
into four modules of about 10x10x40 m'^. An estimate of the costs of the main items of such 
a detector is given in table |^, based on informal contacts with manufacturers. The containers 
and the infrastructure are not included. 



Table 2: Estimate of the cost of the detector 



Item 


Unit cost [CHF] 


Amount 


Price [MCHF] 


Photodetectors 


14.0/channel 


1041667 


14.6 


WLS fibres 


1.0/meter 


10417 km 


10.4 


Liquid Scintillator 


1.0/kg 


15 ktons 


15.0 


Electronics 


10.0/channel 


1041667 


10.4 


Total 


50.4 



3 Event selection and background rejection 

The fiux components of the NGS neutrino beam are reported in table 0. The expected event 
rates for pure quasi-elastic, resonances and deep inelastic processes were computed using the 
latest version of the NGS neutrino beam p5[ (they are given in table ^). The number of Vr 
interactions for several values of Am^ are shown. The rates are given per kiloton assuming four 
years of NGS beam at 4 x 10^^ protons on target per year. 

From table ^, it can be seen that the quasi-elastic processes for oscillated Vj. are a significant 
fraction of the total number of Vr interactions. Furthermore, the resonant interactions, which 
also make an important fraction of the total number of events, can also have a topology similar to 

Among these, the use of an oil tanker is a possibihty. 

^We define the deep-inelastic processes as those which occur above W"^ = 2 GeV^. 
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Table 3: Neutrino beam components at the Gran Sasso location. The numbers indicated with 
a -k are those of the contamination quoted in Reference . The other values are those of the 
latest version of the NGS beam p5[ 



Source 


total flux [z//(pot m^)] 


relative importance 




7.2 X 10-9 


100% 




1.6 X 10-1° 


2.2%^ 




4.4 X 10-" 


0.6% 




5.8 X 10-12 


0.08%* 



the quasi-elastic processes. Therefore, we focus our analysis on the clean quasi-elastic topology 
where the r decays to an electron. We define this topology as an electromagnetic shower and 
at most one additional charged track in the final state. We shall see how it is possible to have 
a good efficiency by a simple selection of these events while rejecting an important fraction of 
the background coming from the u^, contamination of the beam as well as the other sources of 
background. 

The selection of events could be refined by studying other event properties such as missing 
transverse momentum. Crucial experimental data giving information about the properties of 
the background, including cross-section measurements, could in principle come from another 
experiment like 1216, where oscillations are excluded at large mixing angle. 

It should be mentioned that other decay channels of the r can be studied in addition to the 
electronic channel. For instance, the t ^ tt^Ut decay in quasi-elastic interactions was studied 
in CHARM II |2^ for a ^ u^- oscillation search at high sensitivity. This channel could 



provide a signal independent of the u^, contamination of the beam. 

Several processes can produce events that mimic the quasi-elastic topology of the signal. 
The following contaminations were studied: 



• i/e contamination 

The main source of background for an appearance experiment selecting quasi-elastic i^r 
interactions where the r decays to an electron is the contamination of the beam. As 
can be seen in table ^, this contamination is of the order of 0.6% of the z/^ component. 
Figure H(a) shows the spectrum of z/^ and z/g of the NGS beam. 

Figure ^(b) shows the neutrino energy distributions of quasi-elastic events and of 
Ur quasi-elastic events where the r decays to an electron, assuming different values of 
Am^. It can be seen that the neutrino energy distribution of the background from the 
contamination of the beam is harder than the one of the oscillated quasi-elastic events. 

One method to reject this background is to cut on the energy of the identified electron. 
Since the oscillation probability favors the low energy range of the z/^ spectrum, and since 
the electron of the signal comes from the decay of a r, the electrons of the signal will have 
a much lower energy than the electrons of pure quasi-elastic z/e processes. Figure ^(a) 
shows the electron energy distributions from both quasi-elastic z/g interactions and r 
decays from u^- quasi-elastic interactions. Figure |^(b) shows the normalized integral of 
these distributions between zero and a given value of the electron energy (x-axis). A cut at 
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Table 4: Expected number of events per kiloton at 735 km for different processes for four years 
of NGS beam at 4 x 10^^ pots/year. The numbers in parentheses are the numbers of events 
with a T decaying to an electron for the quasi-elastic (QE) and resonant (RES) processes. All 
processes are charged-current interactions except where mentioned (NC) 



Source 


Events/kton 


Vr QE [Am^ = 5.0 X lO^^ eV^; 






Ur QE [Am2 = 2.5 x lO^^ eV^; 




fi 7Q (^ 21 


QE [Am2 = 1.5 X 10-3 eV^' 




4 39 fO 78) 


Ur QE [Am2 = 1.0 X 10"=^ eV^] 


1.11 (0.20) 


Vr RES [Am2 = 5.0 x 10"=^ eV^] 


^3 9 A7\ 


Vr RES [Am^ = 2.5 x lO'^ eV^] 


14 7 l'9 fil "1 


Vr RES [Am2 = 1.5 X 10"=^ eV^] 


9 49 (1 fiQ) 


Vr RES [Am2 = 1.0 X 10-3 eV' 




2 41 (0 A3) 


Ur DIS [Am^ = 5.0 x 10-^ eV^] 


156 


DIS [Am2 = 2.5 x 10-^ eV^] 


40.0 


Vr DIS [Am2 = 1.5 X 10-^ eV^] 


14.5 


Vr DIS [Am2 = 1.0 X 10-3 eV^] 


6.45 


QE+RES 


1060 


RES (NC) 


190 


z/e QE+RES 


6.4 


DIS 


8320 


Z/e DIS 


74 



1 GeV, which will be described later in the text, has been applied before integrating. Since 
the energy distributions of the electron are similar for all the Am^ values considered, the 
four curves of the signal are overlapping and cannot be distinguished. Another important 
feature of this figure is that the integral of the signal rises very rapidly at low energy 
whereas the background rises slowly. A cut on the energy of the electron can thus reduce 
the background efficiently while keeping most of the signal. For instance, for a cut at 
10 GeV, 75% of the signal is kept while 70% of the background is rejected. Figure |^(c) 
shows the ratio of the number of background events and the number of signal events as a 
function of Am^, for a cut at 10 GeV. Figure ^(d) shows the signal selection efficiency for 
a 10 GeV cut as a function of Am^. The final choice of the energy cut should be done by 
maximizing the sensitivity of the experiment, which will critically depend on the choice 
of other event selection criteria, baseline and detector mass. We have checked a posteriori 
that the sensitivity does not change significantly if the maximum allowed lepton energy 
varies in the interval between 4 and 12 GeV (with a 735 km baseline and a 15 kiloton 
detector). Therefore, in the present study, we fix the upper limit of the lepton energy 
to be 10 GeV independently of the baseline. We conclude from figure ^ that a good 
rejection of quasi-elastic z/g events can be achieved by a simple energy cut. Additional 
criteria to reject the contamination could be based on the kinematics of the events (for 
instance, the angle of the lepton with respect to the beam), but have not been used in 
this analysis. 
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• in NC processes 

The conversion of the photons coming from a 7r° decay in a neutral-current process can fake 
the electron signature. To reject this background, we apply a combination of topological 
cuts on the signal recorded in the scintillator. Assuming that the vertex position is known, 
a 7r° can mimic an electron in the following conditions: 

— at least one photon converts in a strip adjacent to the vertex position, while the two 
electromagnetic showers of the photons overlap. 

— one photon converts in a strip adjacent to the vertex position, while the other is 
missed (for example in an asymmetric decay of the 7r°). 

In the cases stated above, the pulse height of the signal recorded in the scintillator can 
be used to discriminate between the passage of a single electron (which will behave as 
a minimum ionizing particle in the first tenths of a radiation length) and two electrons 
coming from a photon conversion. With a statistics of 10 photoelectrons per minimum 
ionizing particle passing through a strip, it is possible to reject 92% of the two-electron 
background using the information of the first strip only (see figure |^). The selection 
efficiency of the signal is then about 86%. 



Given the previous considerations, the selection of Ur events with a quasi-elastic topology 
where the r decays to an electron is the following: 



1. The event must have a pure quasi-elastic topology, with a visible vertex identified by an 
electromagnetic shower and at most one charged track (assumed to be the proton). In 
the case where the proton track is not seen, we require a large pulse height corresponding 
to a heavily ionizing stopping particle. This peak will then be identified as the vertex. 
To ensure a well defined topology, we require that the proton track should not interact. 

2. The electromagnetic shower should be connected to the vertex. 

3. The visible energy of the electromagnetic shower should be between 1 GeV and 10 GeV 
to ensure that it comes from a r decay, and not from a quasi-elastic event. The lower 
cut on the energy also cuts 7r° events, which peak at low energy. 

4. The energy deposit of the electromagnetic shower in the first plane should be compatible 
with a single minimum ionizing particle. 



The selection efficiencies for signal and background were evaluated using a quasi-elastic and 
resonant interaction generator and a Monte Carlo parametrized to simulate the detector 
response and reconstruction. Table ^ shows the different rates of events as computed using 
the cuts mentioned above. The energy distribution of the selected events is shown in figures ^ 
and 1^ for the two baselines under study. Additional criteria could still be explored, as mentioned 
earlier. As an example, figure ^ shows the angular distribution of the lepton with respect to 
the beam direction for the signal and the background. 
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Table 5: Expected number of events for different processes at 735 km for four years of NGS 
beam running at 4 x 10^^ pots/year assuming the efficiencies discussed in the text 



Source 


Events/kton 


efficiency 


Vr QE+RES [Am2 = 5.0 x IQ-^ eV^] 


3.91 


28% 


QE+RES [Am^ = 2.5 x IQ-^ eV^] 


1.07 


28% 


Vr QE+RES [Am2 = 1.5 x IQ-^ eV^] 


0.69 


28% 


Vr QE+RES [Am^ = 1.0 x lO^^ eV^] 


0.18 


28% 


i^e QE+RES 


0.44 


6.9% 


RES (NC) 


0.042 


0.022% 



4 Sensitivity to neutrino oscillations 



The sensitivity is the "average upper hmit that would be obtained by an ensemble of experi- 
ments with the expected background and no true signal" , and has been computed by using the 



statistical techniques reported in reference |28|. The sensitivity of the proposed experiment has 



been computed for a detector of 15 kilotons and a baseline of 735 km. The resulting sensitivity 



plot is shown in figure 10. The systematic uncertainty on the contamination of the beam is 



assumed to be of the order of 5%. 

We studied the possible improvements that can be achieved by varying the mass and the 
baseline. Figure |^ shows the lowest Am^ reachable by the experiment at 735 km as a function 
of its mass. Since the experiment is not background- free, the minimum Am^ value does not 
scale as l/y^niass. 

When Am^(eV^) <^ £'(GeV)/L(km), the oscillation probability V is such that the gain in 
the number of oscillated v-r compensates exactly for the loss of neutrino flux 0(z/^) as the baseline 
is increased: V^u^ z/t-) x0(z/^) ^ constant. On the other hand, the background is proportional 
to the V flux, which decreases like 1/LF'. The minimum Am^ value that can be probed with a 
given experiment having a non-negligible background thus decreases by increasing the distance 



from the neutrino source, at the price of a reduced maximum sensitivity. Figure |TT] shows the 
sensitivity plot for the same detector in the same beam, when it is positioned 3000 km away from 
the neutrino source. It should be stressed that the systematic uncertainties on the background 
(dominated by interactions) are in this case less signiflcant with respect to the 735 km baseline 
experiment, since the signal to background ratio improves by a factor (3000/735)^ ~ 17. 



5 Conclusion 



We believe that the liquid scintillator detector technology could be used to build a large detector 
to search for oscillations in the atmospheric region with a long baseline beam. When- 

ever Am^(eV^) ^ i?(GeV)/L(km), with the current beam designs, the quasi-elastic regime 
has a signiflcant contribution to the overall interaction rate of the oscillated Vr. In this study, 
we have concentrated on a simple selection of quasi-elastic interactions where the r subse- 
quently decays to an electron. A 15 kiloton detector running for four years at 735 km from 
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the neutrino source of a beam similar to the proposed NGS, can probably achieve a minimum 
Am^ of about 1.5x10"^ eV^ at 90% CL. in appearance mode. This result could improve down 
to about 1.1x10"^ eV^ if the baseline were increased to 3000 km. Additional channels in the 
quasi-elastic regime, as well as the analysis of deep inelastic events, could improve the potential 
of the experiment in both appearance and disappearance modes. 
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Energy [GeV] 



Figure 3: Neutrino flux of tfie NGS beam at tlie far location (a). Neutrino energy distribu- 
tion (b) of z/e quasi-elastic events and Vr quasi-elastic events where the r decays to an electron, 
assuming different values of Am^. The normalization of figure b is the same as in table ^ 
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Figure 4: Energy distribution of the electron in quasi-elastic events and in Ur quasi-elastic 
events where the r decays to an electron (a), integral of the distributions of (a) as a function 
of the energy (b), ratio of background to signal for a cut at 10 GeV as a function of Am^ (c), 
efficiency of the 10 GeV cut selection as a function of Am^ 
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Number of photoelectrons 



Figure 5: Probability distribution of the number of photoelectrons assuming an average of 
10 photoelectrons for a minimum ionizing particle crossing a single strip. The peak at lower 
values is for a single particle while the peak at higher values is for two minimum ionizing 
particles 
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Figure 6: Visible energy for the signal and the background at 735 km 
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Figure 7: Visible energy for the signal and the background at 3000 km 




Figure 8: Angular distribution of the lepton with respect to the beam direction for the signal 
and the background at 735 km 
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Figure 9: Minimum Am^ at 90% C.L. as a function of the mass for 735 km 




Figure 10: Sensitivity plot at 735 km 
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Figure 11: Sensitivity plot at 3000 km 
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